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We present the details of a mean-field approximation scheme for the quantum 
mechanics of N DO-branes at finite temperature. The approximation can be 
applied at strong 't Hooft coupling. We find that the resulting entropy is in 
good agreement with the Bekenstein-Hawking entropy of a ten-dimensional 
non-extremal black hole with 0-brane charge. This result is in accord with 
the duality conjectured by Itzhaki, Maldacena, Sonnenschein and Yankielow- 
icz. We study the spectrum of single-string excitations within the quantum 
mechanics, and find evidence for a clear separation between light and heavy 
degrees of freedom. We also present a way of identifying the black hole 
horizon. 



1 Introduction 



The physics of black holes has played a prominent role in our quest to un- 
derstand quantum gravity. Semiclassical considerations have shown that the 
horizon of a black hole has an associated thermodynamic entropy [|TJ , and a 
key test of any proposed theory of quantum gravity should be to provide a 
microscopic explanation of this entropy. 

Dramatic progress was made a few years ago, when certain extremal black 
holes were realized as collections of D-branes in string theory. This descrip- 
tion led to a precise counting of microstates, which was in exact agreement 
with semiclassical black hole thermodynamics @. Unfortunately this count- 
ing relied on supersymmetric non-renormalization theorems, and therefore 
could only be applied to certain classes of extremal black holes. 

A more general understanding of black hole entropy requires a non- 
perturbative definition of string theory. This is now available, at least in 
certain backgrounds, thanks to the M(atrix) and Maldacena conjectures || |J 
(for reviews see ||). These conjectures relate non-perturbative string the- 
ories to dual strongly-coupled large- N gauge theories. In this framework, 
black hole entropy is identified with the entropy of the density matrix which 
describes the gauge theory at finite temperature. 

In principle, one can use these dualities to understand black hole physics 
in terms of gauge theory dynamics. In practice, however, this requires two 
things: a precise map between gravity and gauge theory quantities, and a 
tractable calculational scheme in the gauge theory. Some progress has been 
made on the first issue [0 , although even such basic properties as spacetime 
locality are still obscure from the gauge theory point of view. 

In this paper we will focus on the second issue, of developing practical 
methods for doing gauge theory calculations. The gauge theory is strongly 
coupled whenever semiclassical gravity is valid, so we must study the gauge 
theory non-perturbatively.[] We do this using techniques from self-consistent 
mean field theory. This provides us with an approximation to the density 
matrix which describes the gauge theory at finite temperature. A key test of 
our approximation is whether it reproduces the semiclassical thermodynamics 
of the black hole. As we will see, according to this criterion our approximation 
works quite well, at least over a certain range of temperatures. 

In M(atrix) theory, one can argue that the entropy of certain black holes is not renor- 
malized beyond one loop 0. The argument assumes Lorentz invariance is recovered in 
the large N limit, a property which has been checked at leading order in J8|. 
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For simplicity we will concentrate on the quantum mechanics of iV D0- 
branes, with sixteen supercharges and gauge group SU(N) 0. At large 
N and finite temperature, the effective 't Hooft coupling of the quantum 
mechanics is 
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Note that the quantum mechanics is strongly coupled at low temperature. 
This quantum mechanics is dual to a ten-dimensional non-extremal black 
hole in type IIA supergravity, with N units of 0-brane charge 
metric of the black hole is 
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where c = 2 7 7r 9 / 2 r(7/2) and gyM is the Yang-Mills coupling constant. The 
horizon of the black hole is at U — Uq, which corresponds to a Hawking 
temperature 



T = {9 ™ Nr " 2 <Sf = a2034 (g ™ N) ~ y2 (Sf 
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The dual quantum mechanics is to be taken at the same finite temperature. 
The black hole has a free energy, which arises from its Bekenstein-Hawking 
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Duality predicts that the quantum mechanics should have the same free 
energy. The supergravity description is expected to be valid when the cur- 
vature and the dilaton are small near the black hole horizon. This regime 
corresponds to the 't Hooft large-A limit of the quantum mechanics, when 
the temperature is such that the dimensionless effective coupling (P lies in 
the range [K| 

1 « 91 « iV 10/7 • (5) 



An outline of this paper is as follows. In section 2 we develop a mean- 
field approximation scheme for 0-brane quantum mechanics, building on our 
earlier work |Tj|. In section 3 we present numerical results for the behavior 
of the gauge theory, focusing on thermodynamic quantities. We compare 
our results to the black hole predictions, and find good agreement over a 
certain range of temperatures. Section 4 is devoted to a spectral analysis of 
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the propagators, to extract the spectrum of stretched strings that make up 
the supergravity background. In section 5 we discuss how local spacetime 
physics, such as the size of the black hole horizon, may be extracted from 
the gauge theory. Section 6 gives our conclusions and a discussion of possible 



future directions. A summary of our results has appeared in JT3 



2 Mean-field approximation for O-brane quan- 
tum mechanics 

The basic idea of our approximation is to treat the 0(N 2 ) degrees of free- 
dom appearing in O-brane quantum mechanics as statistically independent, 
with interactions taken into account via a sort of mean-field approximation. 
In the rest of this section we present several reasons to believe this simple 
approximation captures some of the essential physics of the quantum me- 
chanics in the supergravity regime. In the next section we will show that 
the approximation gives results which are in good agreement with black hole 
thermodynamics over a certain temperature range. 

Let us begin by stating our approach to studying strongly-coupled systems 
in rather general terms. We are presented with a strongly-coupled action S, 
in our case the action for O-brane quantum mechanics. We approximate this 
action with a simpler trial action So- All quantities of interest can then be 
computed as an expansion in powers of S — Sq. For instance, the free energy 
has an expansion |L4] 



(3F = (3F - (e-^ 5 °) - l) Cj0 (6) 

= pF + (S-S ) -^((S-S ) 2 )c,o + --- 

where a subscript C, denotes a connected expectation value calculated using 
the trial action So- If the trial action comes sufficiently close to capturing 
the dynamics of the full action, then this expansion should be well-behaved, 
even if the full action S is strongly coupled. 

This sort of approximation relies crucially on an appropriate choice of 
trial action. In our case, we shall take Sq to be the most general quadratic 
action that one can write in terms of the fundamental gauge theory degrees 
of freedom. This means that our trial action involves an infinite number 
of adjustable parameters, namely the momentum-dependent two-point func- 
tions of all the fundamental fields. One can regard these propagators as 
providing an infinite set of variational parameters. To fix these parameters 
we solve a truncated set of Schwinger-Dyson equations. These gap equations 
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provide a non-perturbative approximation to the true two-point functions of 
the theory, by resumming an infinite set of Feynman diagrams. 

As we shall see, this sort of approximation has several attractive features, 
which initially motivated us to apply these techniques to 0-brane quantum 
mechanics. 



• The approximation is non-perturbative in the Yang-Mills coupling con- 
stant, and self-consistently cures the infrared divergences which are 
present in conventional finite-temperature perturbation theory. This 
makes it possible to apply the approximation at strong coupling, at 
temperatures where one can make a direct comparison with black hole 
predictions. 

• We can formulate the approximation in a way which respects 't Hooft 
large- N counting, by only keeping planar contributions to the Schwinger- 
Dyson equations. This means that an overall factor of N 2 in the free 
energy, as well as the appearance of the gauge coupling only in the com- 
bination gy M N, is guaranteed. But this is exactly the form (Q) of the 
supergravity result! That is, we are proposing that the overall factor of 
N 2 in the supergravity free energy can be understood in terms of 0(N 2 ) 
elementary quasiparticles, which are in one-to-one correspondence with 
the degrees of freedom appearing in the fundamental Lagrangian.0 Inci- 
dentally, this means that our approximations are hopeless at couplings 
(outside the range (§)) which are so strong that 't Hooft scaling breaks 
down. 

• A key feature of the approximation is that a quadratic trial action will 
automatically respect all symmetries that act linearly on the fundamen- 
tal fields. This is crucial in a problem like 0-brane quantum mechanics, 
where symmetries play such an important role. By working in a su- 
perfield formalism with off-shell supersymmetry, our trial action will 
have Af = 2 supersymmetry and SO (2) x SO (J) rotational symmetry 
(out of the underlying Af = 16 supersymmetry and 5*0(9) rotational 
symmetry). 

• Another feature of the approximation is that it avoids certain infrared 
problems which are present in the full 0-brane quantum mechanics. The 
difficulty is that the partition function of the full quantum mechanics 
contains an infrared divergence from the regions in moduli space where 
the O-branes are far apart. This leads to a divergent contribution to 
the entropy with an overall coefficient O(N). From the supergravity 

2 A "Fermi liquid" approach to black hole physics! 
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point of view, this corresponds to a thermal gas of gravitons. This di- 
vergence may be regulated by putting the system in a finite box. The 
black hole entropy which is 0(N 2 ) can then easily be made to dominate 
over the O(N) contribution. Our mean-field approximation automati- 
cally computes the 0(N 2 ) piece, while discarding the subleading O(N) 
divergence, so no additional infrared regularization is required. 

This sort of approximation also has some potential drawbacks. 

• An unfortunate fact is that there is no a priori guarantee that the ap- 
proximation works well. One has to choose a trial action and a set of 
gap equations, and hope that with appropriate choices the approxima- 
tion works well. In our case, we will be able to justify our choices a 
posteriori by showing that we get good agreement with black hole ther- 
modynamics over a certain temperature range. Another way to justify 
the approximation is to compute higher-order terms in the expansion 
(||) and show that they are small. We have not attempted this for the 
full 0-brane problem, although toy models show promising behavior 

& 

• Although the approximation respects all symmetries which act linearly 
on the fields, it breaks symmetries that act non-linearly. As there is 
no superspace formulation of theories with 16 supercharges, we can 
only realize a subgroup of the supersymmetries (in our case M = 2) as 
acting linearly on the fields. This is sufficient, for example, to make 
our approximation to the vacuum energy vanish as (3 — > oo. However, 
the remaining supersymmetries and R-symmetries are broken by the 
approximation. Another important symmetry which acts non-linearly 
on fields, and is therefore broken by our approximation, is gauge invari- 
ance. More precisely, our quadratic trial action is not invariant under 
BRST transformations. As we shall show in section 2.3, this difficulty 
can be largely overcome by an appropriate gauge choice. 



Before presenting the details of the approximation, let us note that the 
techniques we are using have a long history. They are closely related to vari- 
ational methods 0] and self-consistent Hartree-Fock approximations, and 
also go by the name of modified perturbation theory [plj. They are equiv- 
alent to the effective action formalism developed in fll6"f . Similar techniques 
have been applied to QCD |T7[ , and related techniques are used to study 
finite-temperature field theory [|l8j. Our own work on the subject began 
with |TIJ, where we were motivated by the 0-brane problem to apply these 
techniques to several toy problems in supersymmetric quantum mechanics. 
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Related techniques have been applied to + dimensional Yang-Mills inte- 
grals in , and have also been used to study Wilson loops in Af = 4 gauge 
theory in p0| . 



2.1 The 0-brane action in M = 2 superspace 

We begin by formulating the 0-brane action in M = 2 superspace. For more 



details see appendix A of |i"2|. 



M = 2 supersymmetry means that we have an SO (2) R-symmetry, with 
spinor indices a, (3 = 1,2 and vector indices i,j = 1,2. The SO{2)r Dirac 
matrices 7L3 are real, symmetric, and traceless. Given two spinors ip and Xi 
there are two invariants one can make, which we denote by 

4> a Xa and 1p a Xa = 2 e «/#aX/3 • 

N = 2 superspace has coordinates (t,9 a ), where 9 a is a collection of real 
Grassmann variables that transform as a spinor of SO(2)r. The simplest 
representation of supersymmetry is a real scalar superfield 

^ = ct> + i^j a + fe 2 . 

It contains a physical real boson <ft and a physical real fermion ip a , along with 
a real auxiliary field /. To describe gauge theory we introduce a real spinor 
connection on superspace T a , with component expansion 

r a = Xa + A 9 a + X^^pOp + deapOp + 2e a/3 X l3 9 2 . 

The fields X % are physical scalars, while \ a are their superpartners, d is an 
auxiliary boson, Xa are auxiliary fermions, and A is the 0+1 dimensional 
gauge field. 

To write a Lagrangian we introduce a supercovariant derivative 

and its gauge-covariant extension 

V a = D a + T a . (8) 

The 0-brane action is built from a collection of seven adjoint scalar multiplets 
$ a that transform in the 7 of a G 2 C 5*0(9) global symmetry, coupled to a 
U(N) gauge multiplet r a . The action reads 

S S YM = ^- [ dtd^Trf—V^^a^-l-V^aV^a-'-U^a^b^c 

9ym J V 4 2 3 

(9) 
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Here Ti = |7lg{V ' a , V/3} is the field strength constructed from T a , and f a bc 
is a totally antisymmetric (^-invariant tensor, normalized to satisfy^ 

3 

f abc f abd = -s cd . (10) 

Strictly speaking, the relative positions of the N 0-branes are governed by 
an SU(N)/Zn gauge theory, but in the large- N limit we can approximate 
this by U (N) . 

We are interested in the finite temperature properties of the action (H) . 
We work in Euclidean space, setting 

S E = ~iS M , t = it, A 0E = -iA 0M , f E = -if M ■ 

Note that we must Wick-rotate the auxiliary fields, to get a Euclidean action 
that is bounded below. As usual we compactify the Euclidean time direction 
on a circle of circumference (3, which is identified with the inverse temper- 
ature. Bosons are periodic while fermions are antiperiodic; for example we 
write the mode expansions 



VP 

1 

7^ 



reZ+1/2 



2.2 Gauge fixing 

Our approximation is based on resumming an infinite class of Feynman dia- 
grams to obtain an approximation for the two-point functions at strong cou- 
pling. To make this procedure well-defined, we must fix a choice of gauge. 
For reasons we will explain, it is extremely advantageous to work in the gauge 

D a T n = 0. (11) 



The first advantage of this gauge is that, since (|TT|) is a condition on su- 
perfields, our gauge choice preserves manifest supersymmetry. In terms of 
component fields, it sets 

d t A = 0, d = 0, \ a = \d tXa - (12) 

This is a complete gauge fixing, i.e. having made this choice there is no 
residual freedom to make additional gauge transformations. 

3 This corrects a normalization error in (131. 
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A second advantage is that our gauge choice is well-defined at finite tem- 
perature. To see this, note that the zero mode of the gauge field, which we 
denote A 00 , survives as a physical degree of freedom. This is important, be- 
cause at finite temperature the gauge theory acquires an additional dynamical 
degree of freedom, namely the value of a Wilson-Polyakov loop around the 
Euclidean time direction U = p e l § dl ~ A a_ \ n our gauge, this physical degree 
of freedom is parameterized by the zero mode. 

U = e^ A °° (13) 

Note that at finite temperature 

A 00 ~ A 0Q + 27r/v^ (14) 

is a periodic variable. 

Corresponding to our choice of gauge we must introduce a ghost action 
(but no gauge fixing term) 

S gh o S t = 4- I dtd 2 6 Tr (D a CV a C) . 
9ymJ 

For the ghost multiplet we adopt the component expansion 

C = a + (3J a + j9 2 

where a and 7 are complex Grassmann fields and (3 a is a complex boson. At 
finite temperature a and 7 are periodic, while /3 is antiperiodic. 



2.3 Slavnov- Taylor identities 

Mean-field methods usually have a difficult time dealing with gauge symme- 
try. The problem is that the Slavnov- Taylor identities are typically violated 
by the approximation. After gauge fixing, Slavnov- Taylor identities arise 
from BRST invariance of the gauge fixed action. BRST transformations act 
non-linearly on fields, but the sort of mean-field approximation we wish to 
use is based on a trial action that is quadratic in the fundamental fields. 
Such a trial action cannot respect a symmetry that acts non-linearly. Thus 
mean-field techniques typically break BRST invariance, and hence violate 
Slavnov- Taylor identities. 

A major advantage of our gauge choice (|TT|) is that many of the Slavnov- 
Taylor identities become trivially satisfied, so that even a quadratic trial 
action can respect many of the consequences of gauge invariance. To illustrate 
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this, we consider a simplified model, which can be obtained from the full 
0-brane quantum mechanics by discarding all fermion and auxiliary fields. 
That is, we study bosonic Yang-Mills quantum mechanics, with the following 
gauge-fixed Euclidean action. 



S = — !— [ drTii-DrXWrX 1 - -[X\X j }[X\X j } 
9ym J ^ ^ 4 

+^(d T A f + d T aD T a} (15) 

Here A is a U(N) gauge field, with D T = d T + i[A , •]. The fields X % are 
adjoint scalars, and a is a ghost field. One subtle point is that the antighost 
zero mode ai=o does not appear in the action, and therefore should not be 
regarded as a true degree of freedom. It is completely decoupled, and any 
correlators involving q.q vanish. 

To illustrate the difficulties with gauge invariance we have adopted a 
general class of gauges parameterized by £. The action is obtained by gauge 
fixing 8 t Aq = f and then functionally integrating over /, with the weight 
exp(— J drf 2 /2g YM £). Our preferred gauge condition d T A = is recovered 



in the limit £ — ► 0. Expanding the fields in Fourier modes, the action (|15[ ) is 
invariant under BRST transformations 



/ 2ttZ 1 . . 

%Aqi = -V \—T a l + ~m i A 0m,a n 
\ P VP m+n=l 



m+n=l 



m+n=l 

= yj a ^ ( 16 ) 

where rj is a Grassmann parameter. Note that the decoupled antighost zero 
mode q.q is indeed invariant under BRST transformations. 

We can use this BRST symmetry to derive Slavnov- Taylor identities in 
the standard way, from the fact that the expectation value of any BRST- 
exact quantity vanishes. For example, we must have 

(6 v (a l A i)} = 0. (17) 

This gives us the following relation among Green's functions. 

(~jAo t -iAoi) + (ai(^j-ai + -±= ^ [A 0m a n ])) = (18) 
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For I = this Slavnov- Taylor identity is trivially satisfied: the first term 
vanishes since I = 0, while the second term vanishes since do is decoupled. 
Here we assume the Aqj two-point function is finite at I — 0. For I ^ the 
second term can be simplified using the following Schwinger-Dyson equation 
(a consequence of the ghost equation of motion) 

(Tr («,«,)> + ~j < Tr &i A o m , On])) = ~9ymN 2 (19) 

^ ' m+n=l 

where for simplicity we have taken a trace to get rid of matrix indices. Thus, 
modulo the use of an equation of motion, the content of the identity (|T7|) is 
the well-known fact that the gauge field propagator at non-zero frequency is 
given exactly by the gauge-fixing term in the classical action: 

(Tr(A ^ ,_ z )) = |^| for/^0. (20) 

In the limit £ —>■ this Slavnov- Taylor identity implies that the modes of 
A with non-zero frequency do not propagate. But this is an automatic 
consequence of adopting our gauge choice (O), which eliminates all non-zero 
modes of A \ In fact, in the full 0-brane quantum mechanics, all Slavnov- 
Taylor identities which just constrain two-point functions are automatically 
satisfied by working in the gauge (\T1\). 

Next let us consider a Slavnov- Taylor identity on a 3-point function. We 
have the requirement 

(5 v (a l X i m X> n )) = 0. (21) 

Using the transformations (|16D , this gives rise to a Slavnov- Taylor identity 
with the schematic form 

-(d T A X i X j ) = (aaX'X 1 ) . (22) 

If the gauge field carries zero frequency this identity turns out to be trivially 
satisfied (for the same reasons that (HD was trivially satisfied at I = 0). If 



the gauge field carries non-zero frequency then this Slavnov- Taylor identity 
is non-trivial. In particular, in the limit £ — > 0, it states that the amplitude 
to emit a gauge boson with non-zero frequency is 0(£). But this property 
is automatically satisfied by working in the gauge flT/[), where the non-zero 
modes of the gauge field are eliminated. Again, the content of the Slavnov- 
Taylor identity (pT|) is automatically taken into account just by working in 
the gauge ([TT]) . 

This pattern is quite general. All non-trivial Slavnov- Taylor identities 
follow from the requirement that correlators of the form 

(*„ («/-••)> (23) 
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vanish. (There must be at least one a, since 5 V increases the ghost number 
by one and you need zero ghost number to have a non- vanishing correlator.) 
If I = this Slavnov- Taylor identity is trivially satisfied. If / 7^ this 
Slavnov- Taylor identity becomes a constraint on a correlators that either 
involve a gauge boson with non-zero frequency, or involve an antighost with 
non- zero frequency. Correlators with A j-^ Q must vanish in the limit £ — > 
0, and this property is guaranteed by working in the gauge ([11]). In fact 
it is not clear to us whether the Slavnov- Taylor identities have any non- 
trivial content in the gauge (|TT]) . In principle it seems that they could give 
constraints on correlators involving antighosts, but at the level of 2-point 
and 3-point functions, no constraints arise which are not already implied by 
the Schwinger-Dyson equations. 

Does this issue of Slavnov- Taylor identities have any practical impor- 
tance? After all, the approximation could work well even though it is not 
gauge invariant. But it turns out that in our case, the gauge choice (|TT|) is 
crucial. We have used mean-field methods to study gauge theories (including 
@, (|15D) in the more general class of gauges, and have found that the 
system of one-loop truncated Schwinger-Dyson equations does not have solu- 
tions when the gauge theory is strongly coupled. We believe this breakdown 
can be related to the fact that the violation of Slavnov- Taylor identities gets 
worse as the coupling increases. 

In any case, at least for 0-brane quantum mechanics, this difficulty can 
be avoided by working in the gauge ([□]). The vertices that appear in the 
gap equations receive no constraints that are not already implied by the 
Schwinger-Dyson equations (quartic vertices that appear in the gap equations 
will not involve a pair of ghosts). The Schwinger-Dyson equations themselves 
will be satisfied at one-loop level, so the approximation is self-consistent. 



2.4 Trial action and gap equations 

In applying mean-field methods to 0-brane quantum mechanics, the first step 
is to choose a trial action. We will adopt the following trial action, which is 
written in terms of component fields expanded in Matsubara modes. 



So = -jTr(U + tf) + J2 ^2 Tr(X;Xl,)-^^Tr( XarXa ,_ r ) 

l a l r r 

+ E ^2 Tr (^>~«) - E 2^ Tr(vc < -) + E ^ Mftr-i) 



J2 -Tricnm) + rl>W - E - Tr (^) (24) 
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Recall that l,m G Z and r,s G Z + | label Fourier modes, a, (3 = 1,2 
are 50(2) R spinor indices, i, j = 1,2 are SO(2) R vector indices, and a,b — 
1, . . . , 7 are indices in the 7 of G 2 . The parameters A, of, . . . can be thought 
of as variational parameters, which we will fix by solving a set of one-loop 
gap equations. 



The action (24) is essentially the most general Gaussian trial action that 
is compatible with the linearly-realized bosonic symmetries of the problem.^ 
Supersymmetry is broken at finite temperature, so we have not imposed 



supersymmetry on the action (|2H|), although as we discuss below supersym- 



metry gets incorporated into our approximation in a natural way. 

There are a few subtle points to note about this action. One point is that, 
due to the periodicity (0), it is not appropriate to adopt a Gaussian trial 
action for A 0Q . Rather we have adopted the unitary one-plaquette model 
action |2~I 



S D = ~Tr(U + U^) 

for the holonomy U = e l ^ A °°. This action undergoes a large- iV phase tran- 
sition when A = 2. As discussed in jTJ], such a transition is expected to 
separate the perturbative gauge theory regime from the super gravity regime, 
presuming couplings to other fields do not turn this into a smooth cross-over. 
A second minor point is that, as discussed in section 2.3, the antighost zero 
mode is not a physical degree of freedom. We have therefore suppressed the 
terms involving a Q in (|24|). 



Corresponding to the action ( |23] ) we have the 2-point correlators 

(^OO^Oo)o = Po (^l^Do — Ol$ l3 $l+m (XarX/3s)o = a r 5 a p5 r+a 

= Aj^V m (Ci)o = g r 5 ab 5 a p5 r+s (ftf b Jo = J5 ab 6 l+rn 

(a/«m)o = sAm ((3 ar (3f3s)o = t r 5 a p6 rs (ia m )o = Ui$lm (25) 



where (• • -) denotes an expectation value computed using S , and where the 
two-point function of the gauge field zero mode is given by 

'H ■ (* A x>2 < 26 > 



pN \ 3 A 



involving a dilogarithm |I2 



Next we need to choose a set of gap equations to fix the parameters that 
appear in our trial action. For most degrees of freedom we will adopt the 



'including the Z2 -R-parity symmetry discussed in appendix A of fll^ 
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one-loop gap equations discussed in |12[ . These equations can be obtained 
by demanding that the quantity^ 

1 2 

hs = fiF + (S n + S IV — S ) — -((S ni ) )c,o (27) 



is stationary with respect to arbitrary variations of the 2-point functions (|25|) , 
where 5*//, Sm, SVv refer to terms in the SYM-plus-ghost action that are 
quadratic, cubic, quartic in the fundamental fields. The explicit expression 
for this effective action is given in appendix A. 

For the gauge field, however, we use a slightly different gap equation. The 
starting point is the Schwinger- Dyson equation for (TrU), which follows from 
demanding that 

(U) = [ dUd{---)Ue~ s (28) 



is invariant under an infinitesimal change of variables U — > gU with g = 
1 + iuj G U(N). At leading order this implies 

<Trt/> = -—(Tr{U^-)) (29) 
VP ° A oo 

where we have dropped higher-order terms in SA 00 , coming from the Campbell- 
Baker-Hausdorff lemma, which do not contribute at one loop in the mean field 
approximation. Evaluating (p9|) to one-loop order gives a relation between 
expectation values computed with respect to S . 

^(TrC/)n = |(-^<Tr(fM o)>n) ■ (30) 

d ( 1 2 

(Siv)o — ■x((Sju) )c,o 



dp 2 V N ' 2 



The relevant one-plaquette expectation values are [21J 
1 f 1 - A/4 A < 2 

iV^ = { V A A>2 (31) 

W WUAa))o _ H(i + i)-Hi-i) 2 Mi-i) a< 2 



N 1-4 A>2 



2A 



while the expressions for (SVv)o and — !((5m) 2 )c,o are gi ven i n appendix A. 
We adopt (|30D as the gap equation that fixes the one-plaquette coupling A. 

5 This quantity can be identified with the two-loop 2PI effective action of Cornwall, 
Jackiw and Tomboulis [Em . 
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Let us pause to note a few important features of this system of gap 
equations. First, in computing expectation values we have kept only pla- 
nar contributions. This means 't Hooft large-iV counting is automatic: the 
free energy will come with an overall factor of N 2 , and the Yang-Mills cou- 
pling will only appear in the combination gy M N. We henceforth adopt units 
which set gy M N = 1, by rescaling all dimensionful quantities as in [Q. Sec- 
ond, since we have consistently worked to one-loop order while including all 
auxiliary fields, these gap equations respect supersymmetry. Of course su- 
persymmetry gets broken at finite temperature, but in the zero-temperature 
limit these symmetry breaking effects go away. 6 Thus as (3 — > oo the bosonic 
and fermionic propagators will be related by supersymmetry Ward identities, 
and the vacuum energy will automatically vanish. 

To summarize, the parameters appearing in our trial action are fixed by 
solving the following set of gap equations. 

i(W)o^(-f(T r( ^)) n )(|E^|E^ 

I r 

+ ^X>?-| B^) v. 2 ) 2 - \ £(^) V) 2 - j £(^) 2 (a?) 2 
-\ £<*) 2 + \ £(^f ) 2 w 2 + ^ £m 2 ) (3 2 ) 

r I r 

1 ,27rL 2 v-^ 9 3z 27rr 14 v-^ .9 2 9 

^ = (-) 2 + ^E^+^E— ^jE^V 

i m 7* in 

-|(^) 2(T i 2 Po - y E 9r9s + ^ E trts ( 33 ) 

^ " " r+s+l=0 " r+s+l=0 



1 ,2irr ^ 3i 2irr v-^ 9 4 2ixr v-^ 27rs 5i 27rr 9 



a r K (3 ' (3(3*-< l f3f3^(3 s 2f3f3 
14i 27rr v-^ . 9 8.27rr. 4 9 14 v-^ , 27rZ. 9A9 



y y i y r+s+ i=o f 



-j E ^ + 7? E E ( 34 ) 

^ r + s +Z=0 ^ r +s+Z=0 r+s+i=0 

1 /27T/.2 4r^ , 4z 2nr 2 2 12 ^ a2 9 

1 ^ r m r r r r r m+n+l=0 



3 Supersymmetry is unbroken at zero temperature in this model. 



14 



12 4.27rZ,, 



r+s+Z=0 



r+ s +l=0 



8 2nl. 

FT- 



(35) 



9 



1 2vrr 12 a2 2 v-^ ,2itL 2a7 

' +- E ^ s ~ ~~b E (^) A ^ 



1 



-- ^ e ^ + 7? E a l9s+g9rP 

^ r+ s +l=0 " r+s+l=0 " 



m+n+i=0 

,2nl 2 1 , 27rZ 2 



r+s+Z=0 



E 2 27r/ q > 

r+s+l=0 



(36) 
(37) 
(38) 



1 _ .2tti I 



1 

Ui 



1 ^ V o 1 ^ — > 1 ^ V ,27?/. 9 1 o 

+ ~R E a ' ts - 2^ E « U ' + —R E ( — ) S <° s + ^A>*r 



/? 



r+s+Z=0 



2/? ^ v (3 



2/3' 



1 \ - 

1 + ? E ^ 



(39) 
(40) 



2.5 Numerical methods 

The gap equations we have described form an infinite set of coupled algebraic 
equations. We now outline the numerical methods that we used to solve these 
equations. For additional details see appendix B of flT2" . 



The first step is to reduce the infinite set of equations (j32|) - (|40|) down 
to a finite set. To do this we use the following asymptotic forms of the 
propagators, which are valid at large momenta. 

° l ~ (2tt///3) 2 + ml ° r ^ 2^7/3 (2vrr//3) 2 + m 2 1 j 

1 „ i2nr/P 2 _ (2tt///3) 2 



(2nl/P) 2 + m 2 A yr ' (2nr/p) 2 + m 2 g 1 (2tt///3) 2 + m\ 
, ~ 1 , „ jgWg „ ^ (2^//3) 2 

' ~ /r> I /mo , o ~ /n /mo , 9 U i 



(2tt///3) 2 + m 2 7 (2vrr//3) 2 + m 2 1 (2tt///3) 2 + m 2 

At leading order these are simply the tree-level propagators.[] Demanding 
7 This reflects the fact that the quantum mechanics is free in the ultraviolet. 
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that these propagators satisfy the gap equations to the first subleading order 
of an expansion in 1/ (momentum) fixes the asymptotic masses to be 



ml = 


9 H Ottt 1 A 

T3Y' + T3^— a ' + l 


E A ?- 


D 2 


m l = 


3 „ 4? 2nr 14 




11 o 

L1 2 

2//° 


m\ = 


4 t — ^ o 12 t — ^ o 6 o 






m 2 g = 


4 % — ^ o 12 % — ^ o 2 o 

^I>? + jE A ?-^ 






ml = 








m 2 s = 


z v 2irr 2 2 






m 2 t = 


1 V - "* 2 1 2 






m l = 


i v-^ 27rr 







The next step is to fix a mode cutoff N. For modes with — N < l,r < N 
we regard the Fourier modes of the propagators themselves as the unknowns, 
while for modes outside this range we parameterize the propagators in terms 
of the eight unknown asymptotic masses appearing in (|4T). The propagators 



with —N < l,r < N are to be found by directly solving the relevant gap 
equations fl3"2]) - (f4"0|), while the asymptotic masses are to be determined by 



solving the system of equations fll^). Note that all these equations are cou- 
pled. For example, we evaluate the high-momentum parts of the loop sums 
that appear in (|32|) - ( |40"D and ([42]) analytically, in terms of the asymptotic 
masses. 

This leaves us with a finite set of equations.^ Our basic strategy is to 
start at high temperatures (3 <C 1, where we have the following approximate 
solution to the gap equations. 

A « 0.418/3 3/2 <rjj« 0.209/3 1/2 A 2 « 0.282/? 1/2 
e 2 « 0.677 M « 1 (43) 



8 8./V + 13 equations, to be precise, where we've taken advantage of the fact that time- 
reversal invariance makes the (bosonic, fermionic) propagators (even, odd) functions of 
their momenta. 
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(All non-zero modes are approximately given by their free- field values.) Then 
we use the Newton-Raphson method p2| to solve the system of equations at 
a sequence of successively lower temperatures. Our numerical results were 
obtained starting at (3 = 0.1, with 

(3 -> min(1.2/5, (3 + 0.25) 

on each step, and with a mode cutoff N = max(3, 5/3). 

Finally, after solving all the gap equations, we wish to compute the free 
energy. This has the expansion given in (El), which we truncate to 

(3F « (3F + (S- S ) - \((S in ) 2 )c fi ■ (44) 

That is, our approximation to (3F is simply the effective action J e g of ( p7|) .P| 
The explicit expression for I e s is given in appendix A. To calculate 7 e fr nu- 
merically we must make use of the asymptotic forms (p|). For example we 
define the following renormalized sum: 

N 

-5>ga? = -E lo gK(( 2 ^//5) 2 + ^)] 

l -N 

+21og(2sinh(/3m (T /2)). (45) 

3 Mean-field results for thermodynamic quan- 
tities 

In principle the trial action we have constructed contains a great deal of 
information about correlation functions in the quantum mechanics. But in 
this section we will just present numerical results for the behavior of three 
basic quantities: the free energy, the Wilson loop, and the mean size of the 
state. 

3.1 Free energy 

At high temperatures, where the gauge theory is weakly coupled, we find 
that the free energy of the system is 

(3F = 6 log/3 + 0(1). (46) 

9 The 0-brane action also has 6-point couplings, but it turns out that (Svi)a — so 
these terms do not contribute to our approximation for the free energy. 
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Figure 1: The solid curve is the power law fit ( Ww for (3F. The data points 
are calculated from numerical solutions to the gap equations. 



This result can be obtained analytically: the gap equations are dominated 
by the bosonic zero modes, and the free energy is dominated by (3Fq. 

In general, for a weakly-coupled theory in + 1 dimensions, one would 
expect the free energy to behave like log/3. But note that, even though 
the gauge theory is weakly coupled at high temperature, the perturbation 
series is afflicted with IR divergences. Thus, to determine the coefficient of 
the logarithm (which depends on the value of the dynamically generated IR 
cutoff) one must re-sum part of the perturbation series. This is a well-known 
phenomenon in finite temperature field theory WE\. In any case, we expect 



a priori that mean-field methods give good results in the high temperature 
regime. 

As the temperature is lowered the behavior of the free energy changes: at 
(3 ~ 0.7 we find that it begins to roll over and fall off as a non-trivial power 
of the temperature. In the range 1 < (3 < 4 the numerical results for the free 
energy are well fit by 

(3F « -0.79 - 2.0 /T L7 . (47) 

This fit to the numerical results is illustrated in Fig. 1. Note that super- 
symmetry is crucial in making such power-law behavior possible. Without 
supersymmetry the free energy would behave as f3F k, f3E in the low tem- 
perature regime ((3 > 1), where Eq is the ground state energy of the system. 

We obtained ( f47|) by performing a Levenberg-Marquardt nonlinear least- 



squares fit |22| to 75 numerical calculations of the free energy, carried out in 
the temperature range 1 < (3 < 4. To estimate the uncertainty in the best 
fit parameters we varied the window of (3 over which the fit was performed 
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(fitting over the ranges 2 < (3 < 4 and 1 < (3 < 3), which leads to: —0.79 ± 
0.06, -2.0 ± 0.1 and -1.7 ± 0.2. 



It is quite remarkable that the power law (fi?D is in excellent agreement 
with the semiclassical black hole prediction pTTf] 

/3F = -4.12 /T L8 °. (48) 

The exponents differ by 6% while the coefficients of the power-law differ by a 
factor of 2. (An additive constant appears in the mean- field approximation 
for the free energy. We will generally ignore this 'ground state degeneracy', 
since it seems to be an artifact of the approximation when applied to systems 
with a continuous spectrum. Similar behavior was noted in ||12|| .) In a toy 



model studied in it was noted that higher order terms in the expansion 
of the free energy @ appear with approximately the same power law depen- 
dence on temperature as the leading term. Thus by computing higher-order 
corrections one might hope for better agreement of the overall coefficient, 
with the power law essentially unchanged. 

As we go to still lower temperatures, we find that the energy d(f3F)/d(3 
calculated in the mean-field approximation begins to drop below the energy 
of the black hole. In fact the mean-field energy becomes negative around 
(3 = 5.8. Ultimately, as f3 — > oo, the mean-field energy does asymptote to 
zero, as required by the M = 2 supersymmetry which is manifest in the 
approximation. But a negative energy clearly reflects some problem with the 
approximation. 

Fortunately, we can be rather precise about exactly where the approxima- 
tion is going wrong: the difficulty is with the Schwinger-Dyson gap equation 
we have been using to fix the value of the one-plaquette coupling A. Although 
we do not know how to write down a better gap equation for A, we can give 
a prescription for fixing A, that will allow us to obtain reasonable results at 
much lower values of the temperature. This may be regarded either as a 
check on our understanding of why the approximation is breaking down, or 
as a way of building a model for the black hole that can be used at lower 
temperatures. Our prescription for fixing A is simply that, when (3 > 2.5 
(the midpoint of our range 1 < (3 < 4), we choose A so that the free energy is 
given by fl4"Tp. The energy E = d(j3F)/d(3 calculated with this prescription 
is shown in Figs. 2 and 3. 



3.2 Wilson loop 

In our approximation the expectation value of the timelike Wilson loop U = 
p e i§drA j g con trolled by the one-plaquette coupling A, as in (^T|). A key 
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Figure 2: Energy vs. (3. For (3 > 2.5 fixing A by fitting (3F to a power law 
leads to the solid middle line, while the Schwinger- Dyson gap equation for 
lambda leads to the lower dashed line. The upper dot-dashed line is the 
semiclassical energy of the black hole. 
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Figure 3: Energy vs. (3. Same as Fig. 2, but plotted on a log-log scale. 
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Figure 4: The one-plaquette coupling A vs. (3. The Gross- Witten transition 
occurs when A = 2. For (3 < 2.5 we use the Schwinger-Dyson gap equation 
to determine A. For (3 > 2.5 the Schwinger-Dyson gap equation gives the 
dashed line, while fitting /3F to a power law gives the solid line. 



feature of the one-plaquette action is that a large- N phase transition occurs 



at A = 2 pTJ . At this value of the coupling the eigenvalues of U spread out 
around a circle, and become sensitive to the fact that the gauge field is a 
periodic variable. It has been argued that just such a phase transition is 
expected to occur in 0-brane quantum mechanics, as the system moves from 



weak coupling into the supergravity regime [12 



Our mean-field results for A are shown in Fig. 4. We present the results 
for A that are obtained by solving the Schwinger-Dyson equation (|30|), as 
well as the results that are obtained from our prescription of fitting f3F to a 
power law. 

Note that in both cases, A increases monotonically with j3. The Gross- 
Witten phase transition takes place when A = 2; with the prescription of 
fitting (3F to a power law this value is reached at (3 = 7.8. Thus, as expected, 
a phase transition takes place as the system moves into the supergravity 
regime |l!|. By adopting the prescription of fitting (3F to a power law, we 
cannot say anything about the order of the phase transition. But if one takes 
the Schwinger-Dyson result for A seriously, then the Gross- Witten transition 
occurs at (3 = 14.2, and is weakly second order (the second derivative of the 
free energy drops by 0.01 in crossing the transition). 

Our prescription for choosing A by fitting (3F to a power law begins to 
break down around (3 = 14, as we find that A rapidly diverges as f3 approaches 
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14.Q By itself, this is not necessarily a problem: infinite A simply means 
that the Wilson loop is uniformly distributed over U(N). But unfortunately, 
we do not have a good prescription for continuing past this temperature. 
Evidently some of the other gap equations (not just the gap equation for A) 
start to break down at this point. Note that this breakdown does not occur 
until well into the strong coupling regime, as an inverse temperature (3 = 14 
corresponds to an effective gauge coupling g% s = /3 3 ~ 3 x 10 3 . 



3.3 Mean size 

Finally, let us comment on the average 'size' of the state. In our approxi- 
mation the scalar fields X l {r) and a (r) are Gaussian random matrices, and 
their eigenvalues obey a Wigner semi-circle distribution. We can define the 
size of the state in terms of the quantities 

I(Tr(X*(r)) 2 > , 

I(Tr(r(r)) 2 ) . (49) 

The radius of the Wigner semi-circle, given by 2y/E?, is shown in Fig. 5. 
Note that the radius stays fairly constant in the region corresponding to the 
black hole. However, because the superfield formalism we are using does not 
respect the full SO (9) invariance, the radius measured in the scalar multiplet 
directions is not the same as the radius measured in the gauge multiplet 
directions. At (3 = 14 we find 

2-Rscalar = 1-81 2i? gaugc = 0.80 . 

This shows that, as expected, the trial action does not respect the under- 
lying 50(9) invariance. Nonetheless, the trial action may provide a useful 
approximate description of the black hole density matrix in the supergravity 
regime. 

In Fig. 5 we have also plotted the Schwarzschild radius of the black hole 

U /2n = 1.89 /T 2/5 . (50) 

Note that, as the temperature decreases, the Schwarzschild radius becomes 
much smaller than the radius of the eigenvalue distributions. It seems ap- 
propriate to identify the radius of the eigenvalue distributions with the size 
of the region U <C (gfyA^) 1 / 3 in which 10-dimensional supergravity is valid 

m 

10 The Schwinger-Dyson gap equation for A has solutions at all temperatures. 
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Figure 5: Range of eigenvalues (radius of the Wigner semi-circle) vs. j3. The 
upper solid curve is for the scalar fields in the scalar multiplets; the lower 
dotted curve is for the scalar fields in the gauge multiplet. The dashed curve 
is the Schwarzschild radius of the black hole. These results were calculated 
with j3F fit to a power law for (3 > 2.5. 



This brings up a subtle issue. The Higgs fields of the gauge theory are 
expected to correspond to spatial coordinates in the supergravity geometry 
(0). But one is always free to reparameterize the radial coordinate in su- 
pergravity. In (Q) we have implicitly made use of the naive identification 
X = U/2n, where X is a Higgs field and U is the supergravity coordinate 
appearing in (Q). This can be justified at zero temperature, because super- 
symmetry fixes the mass of a BPS stretched string in the gauge theory to 
be given by the tree- level formula mw = X, while in supergravity one has 
m w = U/2n fllOfl . However this particular identification is not appropriate 
at finite temperature. A proposal for relating the two radial coordinates has 



been presented in 123] 



An unambiguous way to fix the relation is to use the fact that mw = {U — 
Uo)/2n in the non-extremal black hole geometry. By computing mw in the 
gauge theory, the mapping between the Higgs field X and the supergravity 
coordinate U can be fixed. However one must first take account of the fact 
that one has a continuous distribution of masses in the quantum mechanics. 
The spacetime geometry will only correspond to the lightest of these states 
as we discuss in section 5. This procedure will be studied further in |24l . 
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Figure 6: The 2-point function of the scalar fields in the scalar multiplets, 
calculated at Matsubara frequencies. The solid curve is a fit to the propagator 
with a twin peak ansatz for the density of states. The upper blue curve is 
for p = 3.78, the lower red curve for (3 = 2.03. 

4 Propagators and spectral weights 



Important information about spacetime geometry is encoded in the spec- 
trum of single-string excitations in the quantum mechanics. To extract this 
information from the Euclidean propagators we introduce a spectral repre- 
sentation for the 2-point functions. By inserting complete sets of states, 
one can show that at finite temperature the analog of the Lehmann spectral 
representation takes the form 

<0M0(O)> = f 7"^-^' < r < ,3 (51) 

where the spectral weight p(u) is defined as the thermal average 

m n>m 

We can interpret p{uj)duj as the effective number of single-string microstates 
with a mass between to and uj + du. We will apply this to the scalar fields in 
the scalar multiplets, setting 



poo 

Af= / dwp{u) 
Jo 



(2tt///?) 2 + cj 2 
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Figure 7: The effective density of states arising from the scalar multiplet. 
We show the result for two different temperatures: the blue peaks are for 
P = 3.78, the red for (3 = 2.03. 



In general, solving the inverse problem to extract the density of states 
from A; 2 is a difficult numerical problem, which we analyze in detail in [24 



However some gross features of the spectral density can be easily seen. Con- 
sider the A 2 propagator which is plotted in Fig. 6. At large frequency the 
behavior of the propagator is controlled by the asymptotic mass m\ fl42|) . 
This mass is of order one in 't Hooft units; for example = 2.8 at /3 = 3.78. 
This indicates that states with a mass of order one are present in the spec- 
trum. In addition, note that a clear enhancement of the propagator at small 
frequency compared to its asymptotic form can be seen in Fig. 6. This 
suggests the presence of light states in the spectrum, with a mass of order 
the temperature. These light states make the dominant contribution to the 
entropy. 

To express this in a more quantitative way, we will make an ansatz for 
the form of p(u). The ansatz will allow us to estimate the spectral density 
without performing a complete analysis of the inverse problem. Our ansatz 
is motivated by the results of f24|, where we considered a 0-brane probe of 



the black hole background. The full analysis of pi shows that the density 



of states consists of two narrow peaks (i.e. narrower than a scale of order 
the temperature), one centered at frequency of order the temperature, the 
other at frequency of order the 't Hooft coupling. Motivated by this result, 
we introduce the following ansatz 

, , U)Z\ UZ2 , , 

(u) — mi) z + Y \uj — m 2 ) 2 + Y 
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where we fix the width 7 = 0.02 (any number much smaller than T will do). 
The parameters in (|52"D are determined by fitting to the A 2 propagator. Note 
that the small value of 7 is motivated by pi| ; we could not hope to extract 
such a small width by fitting the ansatz just at Matsubara frequencies. The 
resulting densities of states are shown in Fig. 7. As can be seen from the plot, 
we obtain a very good fit to the propagators at the Matsubara frequencies. 
The fit determines the parameters Zi and with uncertainties at the level 
of a few percent. 

This result clearly indicates the presence of both the high and low fre- 
quency states mentioned above. It also points to a clear separation between 
the two sets of states - a separation which will play an important role in the 
next section. 



5 Resolving spacetime geometry 

In our mean-field approximation, we have modelled the cloud of 0-branes that 
make up the black hole using Gaussian random matrices. As we discussed in 
section 3.3, the eigenvalues of these matrices have very large quantum fluctu- 
ations. Within the gauge theory, we find that the scale of these fluctuations 
is set by the 't Hooft coupling. 

^ 2 = ^(Tr0(r) 2 ) o ~ (g 2 YM N) 2/3 (53) 

In terms of supergravity, this means that the positions of the 0-branes that 
make up the black hole have very large quantum fluctuations. Indeed they 
fluctuate over roughly the entire region (of size U ~ (gy M N) [fLO] ) in which 
supergravity is valid. One might suspect that these large fluctuations are an 
artifact of our approximation, but it has been argued that the scaling (|53"D is 



an intrinsic feature of 0-brane quantum mechanics |25|. 



This raises a very interesting question. How can we recover local space- 
time physics from the quantum mechanics? In particular, given the large 
fluctuations (|53f), how can we resolve the horizon of the black hole? 



The answer is that local spacetime physics only arises as a low energy ap- 
proximation to the quantum mechanics. To recover local spacetime physics 
from the quantum mechanics we must introduce a resolving time, and inte- 
grate out high frequency degrees of freedom^]. The point is that most of the 
TV 2 degrees of freedom in the quantum mechanics have a very large frequency, 

11 We are grateful to Leonard Susskind and Emil Martinec for discussions on this topic. 
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set by the 't Hooft coupling. 



From the supergravity point of view, this energy scale corresponds to the 
energy of a string that stretches across the entire region in which supergravity 
is valid. A low energy observer within supergravity cannot resolve such high- 
frequency fluctuations. Therefore, to recover local spacetime physics from the 
quantum mechanics, we must first introduce a resolving time e, and integrate 
out all modes with frequencies larger than 1/e. With an appropriate choice 
of resolving time, we should recover the expected result, that the 0-branes 
only fluctuate over a region whose size is set by the horizon of the black hole. 

We begin by discussing a single harmonic oscillator. At finite temperature 
the fluctuation in the oscillator position coordinate is 

(x 2 ) ' 



2wtanh(/?cj/2) ' 

We introduce a resolving time, by smearing the Heisenberg picture operators 
over a Lorentzian time interval e. 

x = r J*e*l*x(t) 

' '7T 



The fluctuations in the smeared operators are suppressed when u > 1/e. 

<* 2 > 



2u;tanh(/3o;/2) 

To take this over into 0-brane quantum mechanics, we use the spectral rep- 
resentation flSTp. The fluctuations in the field are given by 

r°° 1 

<*W 2^(^/2) - (54) 

Following our treatment of the harmonic oscillator, we can introduce a re- 
solving time by setting 

roc —w 2 e 2 /2 

U 2 )= / dujpiuj) — -— — — — - . (55) 



One might worry that this definition of horizon radius is ambiguous, since 
it seems to depend on the choice of resolving time. Fortunately, from Fig. 7, 
we see that the spectral density consists of two well-separated peaks. The 
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Figure 8: The dotted curve is the smeared radius of the Wigner semi-circle. 
The solid curve is the Schwarzschild radius of the black hole, as measured in 
the U coordinate. 

low-frequency peak corresponds to states with an energy of order the tem- 
perature; these states can be thermally excited and should be included in 
the fluctuations which make up the horizon. The high-frequency peak corre- 
sponds to states with an energy of order the 't Hooft coupling, states which 
should be integrated out to see agreement with supergravity. Thus any re- 
solving time which keeps the low-frequency peak and integrates out the high- 
frequency peak is acceptable, and will produce the same horizon radius. 

For reasonable values of e we can easily estimate (0 2 ). Rather than use 
the Gaussian cutoff (|55|), it turns out to be more convenient to define a time- 
averaged size by introducing a factor 1/ cosh(/?a;/2) into the integrand of 



(<P 2 )^ dwp{u) . (56) 
J 2lj smh(pu/2) 

The extra factor has the effect of cutting off the integral at to ~ 1//3. This 
corresponds to a reasonable choice of resolving time, e m j3. Defined in this 
way, our estimate for the time-averaged fluctuations in the 0-brane positions 
is simply given by a Euclidean Green's function, c.f. (|51|). 

(0 2 > = (009/2) 0(0)) 

This is easily calculated as a Fourier transform of our momentum-space prop- 
agators. In Figure 6 we plot the resulting smeared radius of the Wigner 
semicircle 



^ = 2a/ (7^ + 2^2 )/9 
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as a function of (3 (we average over the scalar fields in the gauge and scalar 
multiplets).0 The time-averaged fluctuations in the 0-brane positions go 
down with temperature. This is the expected behavior for the size of these 
black holes. The result for R is in rough agreement with the Schwarzschild 
radius U of the black hole (|50|), which we show in the same plot. 

6 Conclusions and future directions 

In this paper we have developed an approximation scheme for 0-brane quan- 
tum mechanics at strong coupling and finite temperature. We presented an 
ansatz for a trial action which captures some of the behavior of the large- 
N quantum mechanics. The parameters appearing in the trial action are 
chosen according to a set of gap equations which resum an infinite set of pla- 
nar diagrams. The approximation automatically respects 't Hooft large-iV 
counting, and also partially respects the supersymmetries and R-symmetries 
of the quantum mechanics. 

Our main result is that we find good agreement with black hole thermo- 
dynamics over the temperature range 1 < (3 < 4. In addition, we studied the 
behavior of a Wilson loop, and found that as expected a large- iV phase tran- 
sition occurs as the system enters the supergravity regime. We also presented 
results on the mean size of the system, and argued that in the supergravity 
regime this mean size (as measured by quantum fluctuations) exceeds the 
Schwarzschild radius of the dual black hole. We could nonetheless recover 
the Schwarzschild radius from the quantum mechanics, by introducing a suit- 
able resolving time; to make our prescription unambiguous it was important 
that the spectral density showed a clear separation between light and heavy 
degrees of freedom. 

We would like to emphasize that in the temperature range 1 < (3 < 4 
our approximation applies strictly to the gauge theory, and makes no use of 
supergravity information. We presented a prescription for fixing the value of 
the Wilson loop, which allowed us to extend the agreement up to (3 — 14. 
The prescription, however, relies on supergravity inputs. 

Our results are based on several technical developments in the use of 
mean-field methods. Some of these developments were reported in our pre- 
vious work [Q. In the present paper, the main new technical problem we 
faced was the difficulty of treating gauge theories using mean-field methods. 

12 The result for R is dominated by the contribution of the zero-frequency Matsubara 
mode. 
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By working in the gauge ([0]), many of the Slavnov- Taylor identities become 
trivially satisfied. This gauge choice was instrumental in enabling us to find 
a consistent set of gap equations, that could be solved in the strong-coupling 
regime. 

There are two perspectives that one could take on this subject. The 
'supergravity' perspective is that, since gauge theory is better understood 
than quantum gravity, we should try to study supergravity phenomena from 
the gauge theory point of view. The 'field theory' perspective is to regard 
O-brane quantum mechanics as an interesting laboratory for developing and 
testing methods to study field theories at strong coupling. 

Depending on which perspective one adopts, there are several interesting 
possible directions for future work. From the field theory point of view, it 
would be interesting to apply mean-field methods to other models, to better 
understand the range of validity of these techniques. Let us mention one 
possibility. One can apply our techniques to pure M = 2 gauge theory, simply 
by dropping the scalar multiplets. The resulting system of gap equations does 
not have a solution in the low temperature regime (3 > 1. Presumably this 
can be related to the fact that the pure gauge model breaks supersymmetry 
spontaneously ||28|| . It would be interesting to understand this connection in 
more detail. 

It would also be interesting to have additional tests of our approxima- 
tion scheme. The supergravity makes further predictions for the behavior 
of the gauge theory, which could be tested. For example, in [^6| a set of 
predictions were made for the scaling exponents of two-point functions of 
certain operators at zero temperature. These were extracted by computing 
Green's functions in the extremal supergravity background, and taking a 
large time/low frequency limit. It was argued that these predictions follow 
from a generalized conformal symmetry that appears in the 't Hooft limit 

li- 
lt would be interesting to test these predictions against our numerical 
results. Even at finite temperature, one would still expect to recover the 
scaling behavior for frequencies satisfying T < w < (gy^jN) 1 ^ 3 . Unfortu- 
nately, the correlators which are predicted to have a scaling behavior involve 
composite operators whose two-point functions we have not yet computed. 
We hope to study this question further in the future. 

Another interesting direction would be to better understand the duality 
between gravity and gauge theory. In particular, it would be interesting to 
understand better how the supergravity properties of spacetime locality and 
causality emerge from the gauge theory. One might hope to see that the 
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horizon of the black hole is reflected in the dynamics of the gauge theory 



along the lines of | 29fl . Also, as we mentioned at in section 3.3, there is the 
subtle question of which radial coordinate in supergravity corresponds to the 
gauge theory Higgs fields. To address these sorts of issues, it is natural to 
introduce a 0-brane to probe the supergravity background. The probe has 
a dual description in terms of a spontaneously broken gauge theory. In [|4j], 
we use mean-field methods to study this problem. 

Ultimately, one might hope to use mean-field methods to study non- 
equilibrium processes in the gauge theory at strong coupling, perhaps using 
some sort of thermofield formalism For example, it would be extremely 
interesting to study scattering of a graviton wavepacket off a black hole. 
Could one see correlations in the outgoing Hawking radiation? 
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A 2PI effective action 



The two-loop, 2PI effective action is defined by 

I eS = (3F Q + (S n + S IV - S ) - ^((S m ) 2 ) c ,o • (57) 

We adopt units which effectively set gy M N = 1, and suppress the overall 
factor of iV 2 in the free energy. Then the free energy of the trial action is 
given by 

PF = /?F D (A)-^loga 2 + ^loga r 

I r 
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- - ^ log Af + 7 ^ log ^ - - ^ log ef 

l r I 

+j2^gs l -2j2 i ogt r + j2 l °& u i ( 58 ) 

1^0 r I 

where the free energy of the one-plaquette model is 

f-f-ilogf + f A<2 

We also have 

+E((f) 2 -?-i) + E(^)V + i 



+ E((^ + 1 ) + 2 ^( ! ^- 1 )-^ ( "'- 1) 



(60) 



z^o 

where the one-plaquette model contribution is 



We also have the contribution of the 4-point couplings, 



2^2 



, , 2 ^-^ 27rr 27rs 3z v-^ 2vTr 2 1 ^-^ r 

<^>° = E — — aA + ^ E — a ^ + ^ E ^ 

^ r,s ^ ^ ^ Z,r ^ ^ l,m 

14 >r-^ 9 9 lAi v-^ . 9 27rr 7 v-^ . 9 9 

E A? ff ?„ + T E ^— a, + - E 

Z,m Z,r Z 

^ + (62) 

where the two-point function of A 00 is defined in (^Bj). Finally, the contribu- 
tion of the three-point couplings is given by 

4<(w>c, = -|E(^)V)Vo - |e(^)V)^ 

14 /27TZ. 2 a2 14 2 14 2 

^ Z+r+s=0 Z+r+s=0 Z+r+s=0 
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